The listing of claims will replace all prior versions, and listings, of claims in the application: 



Listing of Claims: 

1 . (currently amended) A method of obtaining magnetic resonance signals with signal 
separation for at least two chemical species in a heterogeneous magnetic field using rapid 
gradient echo imaging, comprising the steps of: 

a) obtaining first magnetic resonance signals from pixels in an object having at least two 
chemical species using a first repetition time and a first echo time, 

b) obtaining at least second and third magnetic resonant signals from the pixels using 
second and third echo times, wherein the magnetic resonance signals at time, n=l to N, for 
species, i =1 to M, having real, and imaginary, I, parts is given bv: 

M M 

K = +is'„ = ZiPj^Jn -pjdJ+iZipjdjn +PjcJ 

£=1 £=1 

and a least squares fitting of all signals is given by: 

p = (A^A)-'A^S 

where A is a known matrix for M species , 

c) determining a signal estimate for each species and for each pixel by combining all 
measured signals for the pixel using a linear least squares fitting directly on the signals from each 
pixel to decompose the chemical species, assuming a first value of field heterogeneity (y/o), 

d) calculating a first error to the field heterogeneity, wherein error to the field 
heterogeneity is given by: 

y = (B^B) 'B^i 

where B is a matrix given by: 
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where 
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y = [Mlf Ap^ Ap[ Apl Ap[ ... Apl Api^J ^g^ =27rt^Y.i- Pj^jn - PjcJ 
j=i 

M 

e) repeating step c) using the first value of field heterogeneity and the error from step d), 

f) repeating step d) to calculate a second error to the field heterogeneity, and 

g) updating the value of field heterogeneity and repeating steps c) and d) until an 
acceptable error is calculated. 

2. (original) The method as defined by claim 1 wherein M chemical species are present and 
step b) includes obtaining at least M+1 magnetic resonance signals for each pixel. 

3. (original) The method as defined by claiml wherein fat and water are two chemical 
species and step b) includes obtaining three magnetic resonance signals. 

4. (cancelled) 

5. (currently amended) The A method as dofinod by claim 4 of obtaining magnetic 
resonance signals with signal separation for at least two chemical species in a heterogeneous 
magnetic field using rapid gradient echo imaging, comprising the steps of: 

a) obtaining first magnetic resonance signals from pixels in an object having at least two 
chemical species using a first repetition time and a first echo time, 

b) obtaining at least second and third magnetic resonant signals from the pixels using 

second and third echo times, wherein the magnetic resonance signals at time, n=l to N, for 
species, j =1 to M, having real, R, and imaginary, I, parts is given bv: 

K=K+iK=X\Pj^Jn-pjdjnhiTipjdjn+Pj^Jn) 

J=l ^=1 

and a least squares fitting of all signals is given bv: 
p = (A^A)-^A^S 

where A is a known matrix for M species, 

c) determining a signal estimate for each species and for each pixel bv combining all 

measured signals for the pixel using a linear least squares fitting directlv on the signals from each 
pixel to decompose the chemical species, assuming a first value of field heterogeneity (y/a). 

d) calculating a first error to the field heterogeneity. 
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e) repeating step c) using the first value of field heterogeneity and the error from step d), 

f) repeating step d) to calculate a second error to the field heterogeneity, and 

g) updating the yalue of field heterogeneity and repeating steps c) and d) until an 
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where cf =cos(2W^rJ, = sin(2;rA/^fJ, =2^tS- PL - Pfd^ - PfcJ and 
= 27rt^ [p^ + p^c^ - p^fd^ ) are the matrix elements. 

6. (cancelled) 

7. (currently amended) The method as defined by claim 6 1 wherein step a) includes 
obtaining signals from a single coil. 

8. (currently amended) The method as defined by claim 6 1 wherein step a) includes 
obtaining signals from a plurality of coils and steps b) through g) are performed for signals from 
each coil, and further including the step of: 

h) combining field heterogeneity as determined from signals for each coil. 

9. (original) The method as defined by claim 8 wherein field heterogeneity is determined by 
weighting contributions from each coil. 

10. (original) The method as defined by claim 9 wherein the weighting contribution from 
each coil is a function of the square of the magnitude of the image contributed by that coil. 



11. (original) The method as defined by claim 10 where for each pixel, at position r, the 
combined field map is: 
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where P coils collect P independent images. 

12. (original) The method as defined by claim 1 1 wherein the combined field heterogeneity 
from step h) is smoothed by passing through a low pass filter. 

13. (original) The method as defined by claim 11 wherein M images of each species are 
obtained using signals from each of M coils and the combined field heterogeneity, and then 
combining the M images using a square root of the sum of the squares of the images. 

14. (original) The method as defined by claim 1 wherein step a) includes obtaining signals 
from a plurality of coils and steps b) through g) are performed for signals from each coil, and 
further including the step of: 

h) combining field heterogeneity as determined from signals for each coil. 

15. (original) The method as defined by claim 14 wherein field heterogeneity is determined 
by weighting contributions from each coil. 

16. (original) The method as defined by claim 15 wherein the weighting contribution from 
each coil is a function of the square of the magnitude of the image contributed by that coil. 

17. (original) The method as defined by claim 16 where for each pixel, at position r, the 
combined field map is: 

p=i 

where P coils collect P independent images. 

18. (original) The method as defined by claim 17 wherein the combined field heterogeneity 
from step h) is smoothed by passing through a low pass filter. 
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19. (original) The method as defined by claim 17 wherein M images of each species is 
obtained using signals from each of M coils and the combined field heterogeneity, and then 
combining the M images using a square root of the sum of the squares of the images. 
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